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PROCEDURE OUTLINE 


The purpose of this investigation is to develop an analytical method to study the 
vibration characteristics of piezoelectrically forced quartz plates. The procedure 
is schematically shown in Figure 1, and can be summarized as follows. The three 
dimensional governing equations of piezoelectricity, the constitutive equations and 
the strain-displacement relationships are used in deriving the final equations. For 
this purpose, a state vector consisting of stresses and displacements are chosen and 
the above equations are manipulated to obtain the projection of the derivative of 
the state vector with respect to the thickness coordinate on to the state vector 
itself. The solution to the state vector at any plane is then easily obtained in a 
closed form in terms of the state vector quantities at a reference plane. To 
simplify the analysis, simple thickness mode and plane strain approximations are 
used. 
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The governing equations of piezoelectricity consisting of the equations of motion 
and the charge equations of electrostatics are given by Equations (1) and (2) . The 
quantities a ij > u i and ^i are t ^ ie components of stress, mechanical and electrical 
displacements. The constitutive equations are presented in Equations (3) and (4), 
where Cjikx is the elastic stiffness, and e^i, e-y^, E i and s ij are respectively the 
components of mechanical strain, piezoelectric strain constants, electric field and 
dielectric permitivity. The relationship between mechanical strain and displacement, 
and the relationship between electric field and electric potential are given in 
Equations (5) and (6) respectively. 


EQUATIONS OF MOTION 

a ±i.i = P u i,tt 

CHARGE EQUATION OF ELECTROSTATICS 

Di.i = 0 

CONSTITUTIVE EQUATIONS 

CT ij = c ijkl £ kl ' e kij E k 
D i = e ijk £ jk + s ij E j 

£ ij " 05 < u j,i + u i,j> 

Ei = 

Figure 2 
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The plane xl-x3 is taken to be the plane of the plate, and the x2-direction is 
considered as the thickness coordinate. The simple thickness mode approximation, in 
which the various quantities are just functions of the thickness coordinate, is used 
in the analysis. Also, the system is considered to be under plane strain conditions. 
Invoking the above assumptions, and using the contracted notation given by Equation 
(7), the surviving system of equations are presented in Equations (8) through (10). 
Differentiating the last of the equations (10), using the third of Equation (8) and 
integrating the resulting equation twice, the expression for <j> is obtained 
(Equation 11) , where A and B are constants of integration. A constant field does not 
produce any electric field, hence the constant B in Equation (11) is neglected. 
Substituting Equation (11) in Equation (10), the expressions for the non zero stress 
components are obtained, and are given in Equations 12 and 13. 


CONTRACTED NOTATION 


11 , 22 , 


33, 23 OR 32, 31 OR 13, 12 OR 21 




1, 2, 3,4, 5, 6 


(7) 


INVOKING SIMPLE THICKNESS MODE AND PLANE STRAIN ASSUMPTIONS 


a 6,x2 ” P u l,tt 
°2 , x2 - P u 2,tt 
°2,x2 " 0 


( 8 ) 


£ 2 “ u 2,x2 
£ 6 “ 1/2 u l ( x2 
E 2 - -*, x 2 


(9) 


a 2 ” l/ 2 c 26 u l,x2 + c 22 u 2,x2 + e 22 ^,x2 

a 6 ~ c 66 u l,x2 + c 62 u 2,x2 + e 26 ^,x2 

°2 " I/ 2 e 26 u l,x2 + e 22 u 2,x2 * s 22 ^,x2 


( 10 ) 


<t> - { 1/2 e 26 U 1 +e 22 u 2 !/ s 22 + ^2 +• B (11) 

a 2 - a 26 u lx2 + a 22 u 2>x2 + e 22 A (12) 

a 6 - a 66 U 1 , x2 + a 62 u 2 , x2 + e 26 A (13) 


Figure 3 
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A state vector (V) defined by Equation (14) is chosen. The derivatives of the state 
vector with respect to x2 is obtained from Equations (8) through (13) and the 
resulting expressions are given in Equations (17) through (20). The elements of the 
matrices Bl, B2 and B3 are made up of the material constants and derivatives with 
respect to x2 and time (t). 


STATE VECTOR 


(V) - [ { V L } T 

(V 2 ) T 1 t 

(14) 

{ V ]_ } ■= [ U]^ a 2 

] T ; (v 2 } = [ ct 6 u 2 ] T 

{Vl), x2 = [ Bi 

] {V 2 } - A {bx ) 

(15) 

{V 2 },x2 “ [ B 2 

] (Vi) - A (b 2 ) 

(16) 

{V} ,x2 - [ B3 

] (V) - A tb 3 } 

(17) 



l/ a 66 "( a 62/ a 66) d/dx 2 


e 26/ a 66 

Bl - 

0 p d 2 /dt 2 

bl = 

0 


p d 2 /dt 2 0 


0 

B2 = 

- ( a 26/ a 2 2 ) d/dx 2 l/a 22 

b2 = 

e 22/ a 22 


a 62 = c 62 +e 26 e 22/ s 22 > a 66 “ l/ 2 (C66 + e 26 2 / s 22) > 

a 26 “ l/ 2 ( c 26 + e 22 e 26/ s 22 > J a 22 “ c 22 + e 22 2 / s 22 


( 20 ) 


Figure 4 
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A solution to the differential equation given in Equation (11) can be easily 
obtained and is given in Equation (21) , where {VO} is the state vector evaluated at 
x2=0 . The analyst has the flexibility of choosing any plane as the appropriate 
reference plane. The exponential term in Equation (21) can be expressed in an 
infinite series, and the powers of the matrix B3 can conveniently grouped as shown 
in Equation (23). 


v = e B3 x2 {VO} + B3* 1 A b3 
(VO) = {V} x2=0 


B3x2 


{VO] 


I + B3 


x2 + (B3 x2) 2 /2!+(B3 x2) 3 /3!+ j 


{VO] 


( 22 ) 


B3 2 - 

P 

0 1 

B3 3 = 

0 B1Q i 


- 0 

Q - 


B2P 0 - 

B3 4 = 

r p2 

0 1 

B3 5 = 

r 0 B1Q 2 


- 0 

Q 2 J 


L B2P 2 0 


(23) 


[ P ] - [Bl] [B2] 
[ Q ] = [B2 ] [Bl] 


Figure 5 
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Using the expressions given in Equation (23), the infinite series expansion for the 
exponential term can be conveniently grouped as shown in Equation (24). The elements 
present in Equation (24) can be recognized as a convergent series. The resul an 
expression is given by Equation (25). Substituting this expression in Equation (21), 
the final equation for the state vector at any reference plane m terms of the state 
vector at a reference plane is obtained (Equation 26) . 



0 B1Q 2 
- B2P 2 0 


x2 5 /5! + 


(VO) 


(24) 


cosh (x 27P) B1 yQ sinh (x27Q) 

B27P sinh (x2,/P) cosh (x2/Q) 


{VO} 


(25) 


= [R] {VO} 

{V} = [R] {VO} + A [B3]" 1 (b3 } 


(26) 


Figure 6 
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Symbols and Abbreviations 


A, B Constants of integration 
a ij Constant coefficients 
c ijkl Elastic stiffness 

D-[ Components of electric displacement 

e kij Components of piezoelectric strain constant 

Efc Components of electric field 

Sij Components of dielectric permittivity 

u i Components of mechanical displacement 

c kl Components of strain 

4> Electric potential 

P Mass density 

a i] Stress components 
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